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Abstract. 

The EMC effect is studied in the perturbative QCD hard pomeron approach. In the limit 
x — > and for a nucleus with a constant density the effect is found to be a function of a 
single variable which combines its x-, Q- and A-dependence. At large Q 2 the effect dies out 
as \nQ/Q. At small but finite x a change in the anomalous dimension leads to a weaker 
Q-dependence, which almost disappears as x grows and the anomalous dimension becomes 
small. In this approach the nuclear structure function can be expressed directly in terms 
of the proton one. Calculations of the EMC effect on these lines give results which agree 
reasonably well to the existing experimental data for x < 0.05. 
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1 Introduction 

The nature of the EMC effect at small x is well understood qualitatively. At high enough 
Q 2 the incoming photon splits into a qq pair which interacts with the nucleus target in the 
manner typical for a hadronic projectile and thus experiences absorption which makes the 
structure function (per nucleon) become smaller than for a free nucleon [1-5]. An alternative, 
though equivalent, explanation (appropriate for the reference system in which the nucleus is 
moving fast) is that the gluonic clouds created by different nucleons of the nucleus begin to 
overlap and the gluons recombine, so that their density becomes smaller [6-8]. However, as 
far as quantitative predictions are concerned, published results heavily depend on additional 
semiphenomenological assumptions and often contradict each other. In particular the in- 
dependence of the effect results different in different approaches. In a serie of papers [3-5] 
it is argued that qq configurations of a large dimension give the dominant contribution to 
the absorption, which results essentially independent of Q. On the other hand in the gluon 
recombination approach of [6-8] the absorption is obtained as a clear higher-twist effect dying 
out at large Q. 

In this note we calculate the nuclear structure function in perturbative QCD which seems 
to be applicable to the low x region. To simplify the treatment we additionally take the 
number of colours N c to be very large. In this limit the low x nuclear structure function 
can be calculated in a closed form and results quite simple: it is given by the exchange of 
any number of BFKL pomerons [9] between the projectile qq pair and nucleons of the target. 
Corrections have the relative order of 1/iV 2 , which may be not too bad for the physical case 
N c = 3. 

A peculiar property of the BFKL pomeron is that, due to the unperturbative rise of the 
anomalous dimension up to unity, the effective cross section of the coloured dipole formed 
by the qq pair is proportional to its transverse dimension r rather than to r 2 , as taken in 
[3-5]. As a result, a strong violation of scaling is predicted for the proton structure function 
F2 P (x, Q 2 ) at low x, which should rise linearly with Q. 

For the nucleus our calculations reveal that absorptive corrections due to n interactions 
with the target (n > 2) are indeed independent of Q, as stated in [1], and have the order 
m 2 ~ n in the effective quark mass. However for a heavy nucleus contributions with different 
n obtained in this manner cancel in the sum, so that the resulting absorption results inde- 
pendent of the quark mass. A careful study of this result shows that contrary to the afore 
mentioned assertions large spatial separations of the qq pair are not singled out, its typical 
dimension being small and independent of Q, of the order A^^-^Rn, where Rn is the nucleon 
radius. For a nucleus with a constant density inside a sphere of a radius Ra = A^-^Rq the 
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calculated EMC effect turns out to be a function of a single variable 

— 1/3 ^) a - 

which combines its x-, Q- and A-dependence. The effect goes down as \nQ/Q at high Q 2 
and rises with A and 1/x as expected, its absolute magnitude remaining relatively small up 
to x ~ 1CT 6 and not too low Q 2 . 

The asymptotic pomeron approach may only be applicable at quite small x < 0.001, at 
which the existing experimental data on F2 P (x, Q 2 ) do not contradict the predicted linear rise 
with Q. At larger x the growth with Q is less pronounced, which indicates that the anomalous 
dimension of the pomeron becomes smaller and tends to zero at finite x. It turns out that 
the experimental data on F2 P (x, Q 2 ) up to x < 0.05 can be well described by a subasymptotic 
pomeron with the effective anomalous dimension going down logarithmically in x. With this 
subasymptotic pomeron as an input, we have calculated the nucleon structure functions for 
x < 0.05. Comparison to the existing experimental data shows a quite reasonable agreement. 

As a result of a changing anomalous dimension the Q dependence of the EMC effect 
results also x-dependent, becoming weaker as x grows. With the anomalous dimension close 
to zero, the effect becomes independent of Q, in agreement with the statements in [3-5]. 

2 The 7*^4 scattering in the leading approximation in 1/N C 

As is well known, in the high-colour limit all contributions to the scattering amplitude in 
perturbative QCD can be classified according to the topological structure of the corresponding 
Feynman diagrams [10,11]. The leading diagrams have the topological structure of a cylinder. 
For 7*-hadron scattering they sum into a single BFKL pomeron [12]. Diagrams with more 
complicated topologies correspond to multipomeron exchanges with all sort of interaction 
between the pomerons. 

For 7*-nucleus scattering the leading diagrams in 1/N C correspond to the exchange of 
several pomerons with their possible branchings as they propagate to the nucleus (generalized 
fan diagrams, Fig. 1). In the perturbative domain 

g 2 N c « 1 (1) 

only the simplest of them survive at high energies, those without branchings (Fig. la). Indeed 
the energy dependence of the two diagrams of Figs, la and lb is given by (siS2) 2A and s A 
respectively, so that the contribution from Fig. la clearly dominates if s\ is large. If, on the 
other hand, s\ is finite then it continues to be dominant because in this case the contribution 
of Fig. lb includes an extra factor g 2 N c not compensated by a large Ins. 
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Thus in the dominant approximation in 1/N C we have to sum only the diagrams which 

correspond to Fig. la with any number of pomerons. To do this we have to know the coupling 

of the projectile photon to an arbitrary (even) number of gluons, which form colourless pairs. 

In the high-colour limit and in the lowest order in g 2 N c this vertex has been calculated in 

[12]. For 2n gluons with momenta qi, i = 1, 2n, J2 Qi = 0, it is given by 

/ ?v 2 - i \ n r 1 r 2n 
F {2n) (q t ±) = ig 2n [~^^J J o d» J d 2 r p(a, r) JlJexpttm) ~ !]• (2) 

Here p(a, r) is the (dipole) colour density created by the projectile photon as it splits into 
a qq pair with a part a of its longitudinal momentum carried by the quark and a transverse 
distance r between the quark and the antiquark. For longitudinal (L) and transverse (T) 
photons of virtuality Q 2 this density has been calculated in [3]: 

PL(a, r) = £ Z 2 [a(l - a)] 2 K 2 (e f r), (3) 

e 2 N Nf 

(*{«, r) = j^J2 Z 2 {m 2 K 2 (e f r) + [a 2 + (1 - a) 2 ]e 2 K 2 (e f r)}, (4) 

where the summation goes over flavours, 

e) = a(l - a)Q 2 + m) (5) 

and rrif and Zf are the mass and charge of the quark of flavour /. 

The 2n gluons coupled to the vertex (2) form n pomerons as a result of their interaction. 
We assume that the pomerons are formed by pairs of gluons (1, n + 1), (2, n + 2) and so on. 
Upon reaching the nucleus each pomeron is coupled to a nucleon with a vertex analogous to 
(2) with n = 1, depending on the primed momenta q[ and q' n+i , i = 1, ...,n, and involving 
an unknown (nonperturbative) colour density for the nucleon pN{a' i ,r' i ). It is also clear that 
each pomeron has its total momentum equal to zero: % + q n+ i = q\ + q' n+i = 0. Then its 
propagator will be given by the BFKL Green function for the forward scattering [13]: 

ifiy„qi,^/{Nl -1), i = l,...,n, (6) 

where v = pq, q n+ i = —qi, q' n+i = —q\. The factor l/(N 2 — 1) comes from the projection 
onto the colourless state. 

Integrations over q, L and q[ with the exponential factors contained in the vertices transform 
this Green function to the coordinate space. As a function of coordinates, the BFKL Green 
function vanishes at the origin, so that terms with unity instead of exp(%r) in (2) give no 
contribution. As a result each pomeron propagator turns into the coordinate Green function 
for the forward scattering 

4i/(z/,r,rO/(iV c 2 -l), i = l,...,n, (7) 
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which should be integrated over r and all r\ with the corresponding colour densities. 

One has to finally recall that each interaction with a nucleon of the target is accompanied 
by a nuclear profile function factor T(b) giving the probability to meet a nucleon at a given 
impact parameter b. Also a binomial factor C\ should be included for n interactions. 

Then in the end we arrive at the following expression for the 7* : A forward scattering 
amplitude with n interactions: 

A^(p,q) = -MvC n A w n da J d 2 r p(a, r)[-a(r)/2] n . (8) 

Here 

a{r) = (2^p iVf I da J d2r 'PN(*'y)f(^,r>) (9) 
has the meaning of the cross section for the scattering of a colour dipole formed by the qq 
pair of transverse dimension r off the nucleon. The factor w n gives the total probability to 
find n nucleons at the same impact parameter: 

W n = J d 2 bT n (b). (10) 

Evidently the amplitude (8) has an eikonal form under the sign of integration over a and 
r. The 7M total cross section is obtained as 1771^4/(2^). From (8), 

<7$° = -2Clw n j\a J d 2 r p(a,r)[-a(r)/2] n . (11) 

Summing over all possible numbers of interactions we find the total cross section: 

<T A = 2 J d 2 b J\a J d 2 r p(a,r){l- [1- (l/2)a(r)T(b)] A } . (12) 

Now we turn to the cross section <r{r) defined by (9). The explicit form of the BFKL 
Green function which enters (9) is [13] 



(13) 



where 

oj(k) = -(g 2 N c /27r 2 )[Re ^(1/2 + in) - V(l)]- (14) 

We have retained only the isotropic part, which gives the dominant contribution at high v. 
At such v small values of n contribute in (13) for which lo(k) has the form 

w(k) = A-ok 2 , (15) 

where, in terms of a s = g 2 /4:ir, 

A = (a s iY c /vr)41n2 (16) 
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is the pomeron intercept and 

a = (a s N c /7r)U((3). (17) 
With (15) the asymptotical expression for the Green function results 



/(W'^oo ^ 2rrV A ,/-— exp f- 1 " ^ j . (18) 

As one may see from (8) and (9) characteristic values of r and r' are of the order 1/Q 
and 1/mjv respectively, where m^r is the nucleon mass. Then the exponential factor in (18) 
is close to unity, unless Q becomes so large that In 2 Q ~ In is, which we assume not to be the 
case. Omitting this factor we obtain 

a{r) = prR N , (19) 

where 

R N = J da j d 2 r'r'p N (a',r') (20) 

is the average transverse dimension of the nucleon and the factor jjl is given by 

= l^al 12 N 2 c -l u A 
M v/14iV c C(3) iV 2 VrnT ' 1 ' 

It becomes large as v —>■ oo. 

Eqs. (12) and (19) allow in principle to calculate the 7M cross section in a straightforward 
manner for the longitudinal and transverse photons using the densities given by (3) and (4). 
The nuclear structure function at low x can then be calculated as 



Q 2 



F 2A = ^ [a ( P + af \ . (22) 

However before doing so, we have to make two important remarks. 

The first one concerns the dependence of <r(r) on r. It is proportional to r rather than 
to r 2 as assumed in [3-5]. This fact can be traced to the anomalous behaviour of the BFKL 
Green function at small r, resulting in its turn from the nonperturbative rise of the anomalous 
dimension to unity around the BFKL singularity point. As a result, a large violation of scaling 
is predicted for the structure function at small x: it grows linearly with Q. This fact will be 
clearly seen in what follows. 

Our second remarks concerns the validity of the expression (12). It has been obtained 
under the standard assumption that the radius of the strong interaction is much smaller than 
that of the nucleus. Technically it reduces to the strong inequality 



a(r)T(b) « 1. 



(23) 
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In fact, if cr(r) ~ r%, where tq is the effective interaction radius, then (23) tells us that 
ro << Ra- With (19) and (3)-(5) one observes that a(r) can take on large values for light 
quarks, of the order 1/rrif, when the standard assumption that the projectile can interact with 
only one nucleon at a time becomes wrong. However, on physical grounds the expressions (3) 
and (4) cannot be true for very small values of e/, which allow for very large values of r. At 
such distances nonperturbative confinement effects should become notable, which prevent the 
qq pair from having too large dimensions, larger than some maximal possible dimension v max 
of the order of a typical hadronic size. Then Eq. (23) will always be fulfilled. To simulate the 
confinement effects we shall substitute the small masses mj of light quarks by an effective 
regulator mass m ~ l/r max . 

Finally, to simplify the treatment, for A » 1 and using (23) we can, as usual, transform 
(12) into an equivalent expression: 

o A = lJ d 2 b£da J d 2 r p(a, r) {1 - exp[-(l/2)Aa(r)T(b)]} . (24) 

3 The nuclear structure function in the limit x — > 

We begin with the main part of the structure function given by the cross section for the 
transversal photon. Putting the transversal photon density (4) and the cross section a(r) 
given by (19) into (24) and changing the variable r ^ r/e/ we obtain 

,m _ 2 ^ ? * f a f' ^ f ,2 r 



2 

n ' f r K 2 (r) + [a 2 + (1 - a) 2 ]K 2 (r 



f= 

AfirR N T(b) 



1 — exp 



2 y Ja(l-a)Q 2 + m 2 f/ 



a(l - a)Q 2 + mj 

"(25) 

Expression (25) depends on the quark masses rrif. From its form it is evident that the 

IT) 

dominant contribution to a A comes from the region of integration in a where 

a(l - a)Q 2 < {A^R N /R 2 A ) 2 ~ {u^A^/Rn) 2 . (26) 

In the asymptotic region x — > and/or A » 1, to which this and the following sections 
are devoted, v^A 1 ^ » 1, so that the right-hand side of (26) is large. Then large values 
of a(l — a)Q 2 will give the dominant contribution, which correspond to small transverse 
dimensions of the qq pair of the order 

r ~ Rn/^A 1 / 3 ). (27) 
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Evidently the masses do not play any role at such small values of r, so that we can safely 
put nif = in (25). Then it simplifies to 



= J d% J\a J d 2 r[a 2 + {l-a) 2 ]K 2 {r) 



(2vr 

where 



f AfirR N T(b) 
1 — exp 1 



2Q y /a{l-a) i 



, (28) 



N 



f 



E Z ) ( 29 ) 
/=i 



is the average quark charge squared. 

Note that the cross section on the free nucleon follows from (28) in the lowest order in 
T(b) (and divided by A): 

(T) _ e 2 N C Z 2 yR N yi q2 + (l- a) 2 , 

a ^ - (2vr)3 q y da y^r^) y d rrKi(r) - (30) 

Calculating the integrals over a and r, 

(T) 9vre 2 N c Z 2 fiR N 
° N = 5T2Q • (31) 

It is evident that at high Q the cross section a^p ~ so that the EMC effect dies out. 

It is instructive to compare the total cross section (28) with separate contributions coming 
from a given number of interactions n. From (11) we get: 

N, 



a 



(T) „) 2e 2 iV c 
A (2vr)3 



■C\w n J^Zj [ da [ d 2 r 

j ^ J J 



m 2 



a(l - a)Q 2 + m| 



4— ^K 2 (r) + [a 2 + (l-a) 2 ]K 2 (r) 



(32) 



2y / a(l - a)Q 2 + m} y 

One observes that all terms except the first one diverge as rrif — > 0. The n = 1 term is 
independent of mj in this limit and behaves like l/Q (see (31)). Other terms have the order 

ap 2) ~ Q- 2 ln(Q 2 /m 2 ), " >2) ~ <T 2 m 2 ~". (33) 

Thus they are all of the same order Q~ 2 with respect to Q and therefore asymptotically 
smaller than the dominant term with n = 1. For the structure function (22) they give 
contributions essentially independent of Q. In this respect the conclusions in [5] are correct. 
However comparison to the total cross section (28), which is independent of m/, demonstrates 
that the rescattering terms taken separately have little relation to the physical cross section. 
Evidently all their terms leading in 1/m/ cancel in the sum. The final expression (28) does 
not admit developing in the multiple scattering series. 
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To simplify the calculation of the cross section (28) we choose a simplified nuclear profile 
function T(b) corresponding to a finite nucleus with a constant density 

T(b) = (2/V A )^R 2 A -b 2 , (34) 

where Ra = A 1 /^Rq and Va is the nuclear volume. Then the integration over b can be done 
explicitly and we obtain 

a ( p = (2e 2 N c Z 2 R 2 A /^) J T , (35) 
where Jt is given by an integral over a and r: 

J T = da rdra 2 Kl{r)x{z). (36) 
J o J o 



Here 

1 1 

2 \~z + ~z 



x{z) = - n -- 2 +e- Z [- + - 2 ) (37) 



and 

3 A^jxtRn _ r 

The cross section for the longitudinal photon can be studied in a similar manner with the 
colour density given by (3). In the limit mj — > it is given by a formula analogous to (28): 

AfirR N T(b) ' 



(L) 8e 2 N c Z 2 f l2h f 1 f 2 2 

= (2 ) 3 y J J ' _ a )Ko( r ) 



1 — exp 



2Q v /«(l-«) / 



(39) 



(T) 

and possesses the same properties and behavior in Q as a A . For a free nucleon the longitu- 
dinal cross section turns out to be 

(L) ne 2 N c Z 2 fiR N 

= 256Q • (40) 

Finally, with the nuclear profile function chosen according to (34) we obtain 

4 L) = (4e 2 N c Z 2 R 2 A /7r) J L , (41) 

where now 

pi roo 

J L = da rdra(l - a)K 2 {r) X (z), (42) 
Jo Jo 

with x(z) and z defined by (37) and (38). 

According to (35)-(38) and (41), (42) the EMC effect turns out to depend on a single 
variable zq which combines its x-, Q- and A-dependence. Indeed for the ratio R^^XjQ 2 ) of 
the nuclear structure function per nucleon to the proton one we find 
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where, according to (38), 

3 A^fiRN 

z ° = T,^rJ~ ■ (44) 

In principle, zq may take on arbitrary values, from very small ones up to very large ones. 
With the growth of zo the effect (i.e., 1 — R^) also grows from zero to unity (Fig. 2). At 
small zq it behaves like zq lnzo (see Appendix), which corresponds to going down like InQ/Q 
at high Q 2 . 

To find physically relevant values of zq we note that using (22), (31) and (40) we find the 
proton structure function: 

F 2p (x, Q 2 ) = (ll/512)N c Z 2 fiR N Q. (45) 

As mentioned in the Introduction, it exhibits a strong violation of scaling: it rises linearly 
with Q, which means that the anomalous dimension has risen from its presumably small 
value up to unity in the asymptotical region. Eq. (45) may be used to express the unknown 
product uRn in (44) via F2 P (x,Q 2 ) to relate the parameter zq to the observable proton 
structure function: 

Z0 = T^N c Z 2 Q 2 R 2 F2p{x ' Q2) - 
This allows to calculate the nuclear structure function entirely in terms of the proton one 
without any unknown parameter, 

Here the integrals Jt,l are calculated according to (36) and (42) with zq given by (46). 

To apply these asymptotic formulas to realistic nuclear structure functions at given x and 
Q 2 we have first to check that the proton structure function is reasonably well described by 
the asymptotic expression (45). If we take u ~ 1/x in /j, (Eq. (21)) then it leads to 

F 2p (x,Q 2 ) = c f X , (48) 

with 



384 - 4 " 3 - ■ -< (46) 



c= J — c „ o?J 2 Z 2 R N . 49 

32 V 14C(3) N 2 s K ' 

Comparing (48) with the experimental behaviour of F2 P (x, Q 2 ) one may hopefully determine 

A and c and from (16) and (49) find the parameters a s and Rjy. (Note, however, that the 

exact scale factor which divides v in expressions like i/ A or In v cannot be determined in the 

lowest order BFKL theory. Therefore the linear dependence on Q in (48) may be modified 

by terms of the order a s lnQ, which are assumed small in the theory, but are not so small 

for realistic values of a s .) 
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Studying the existing experimental data for F2 P (x,Q 2 ) at low x [14-17], we find that the 
simple formula (48) describes them quite well for x < 0.001 (Fig. 3), with the parameters A 
and c chosen according to 

A = 0.377, c = 0.0536. (50) 

They correspond to a s = 0.14 and Rn = 0.44 fm. Thus the asymptotic hard pomeron, 
at least, does not contradict the experimental data for the proton structure function for 
x < 0.001 and leads to physically reasonable values of a s and Rn. This provides a justification 
to apply it to the nuclear structure function in the same region. 

Taking (48) for F^p and using (46) and (47) with four flavours and Rq = 1.35 fm we have 
calculated the ratios R^ for 10~ 6 < x < 10~ 3 , Q 2 = 2.5 100 GeV 2 and A = 12^-208. Our 
results are shown by solid curves in Figs. 4-6, which illustrate the x-, Q- and ^4-dependence 
of the EMC effect respectively. The effect rises with 1/x due to the rise of F2 P , falls with Q 
as In Q/Q and, naturally, rises with A as A 1 / 3 . Its overall magnitude remains not very large 
in the whole region explored (1 - < 0.5). 

Unfortunately there are no experimental data on the nuclear structure functions at such 
low x. To be able to compare our results with the existing data, we thus have to move to 
larger x, where the asymptotic formula (48) does not work. 

4 Realistic x and A 

In this section we want to consider the region of not so small x and not so large A, where 
neither the asymptotic form (48) for the proton structure function is valid nor the right-hand 
side of (26) is large. 

The latter condition can easily be taken into account retaining in e 2 , Eq. (5), the quark 
masses squared m 2 (with those for the light quarks substituted by the regulator mass squared 
to 2 ). Now both terms in ej (Eq. (5)) become comparable. It means, of course, that now the 
dominant configurations of the qq pair may reach a size of the order 1/to, i.e., a typically 
hadronic size. 

As to the proton structure function, we have to take into account that its rise with Q 
becomes considerably slower at not so small x, which implies that the effective anomalous 
dimension of the pomeron is, in fact, ^-dependent and goes down with the growth of x. Of 
course one may try to study the behaviour of the pomeron at lower energies (higher x) using 
the BFKL equation with some appropriately chosen boundary conditions (or a driving term) . 
However, since these are poorly known, the resulting predictions are not reliable. Besides, at 
not too small x subdominant terms in 1/lnx may become essential, which are beyond any 
control (see however [18], where such terms are introduced from some theoretical considera- 
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tions). For these reasons, rather than try to determine the subasymptotical behaviour of the 
pomeron on (poorly known) theoretical grounds, we use a simple extension of the asymptotic 
formula (48) to finite x which only takes into account a change in the anomalous dimension, 
to be determined from the comparison to the experimental data. Namely, we assume that 
the effective anomalous dimension of the pomeron is x dependent and smaller than unity at 
finite x. Correspondingly we take instead of (19) 

a(r)=^)(m N r)^/m 2 N , (51) 

where /3(oo) = 1 and (3{v) > 1 at finite v (perturbatively (3 ~ 2). We have also introduced 
the nucleon mass in (51) as a natural scale. Evidently as v — > oo (51) goes into (19) with 

£(z/) = [i(v)m N R N . (52) 

With the cross section u{r) chosen according to (51) and the quark masses rrif retained, 
our formulas for the nuclear cross sections (35)-(38), (41) and (42) slightly change because of 
the flavour dependence. Instead of (35) and (41) we now find 

N f 

a (T,L) = {2 n e 2 NcR 2 AM £ z p(J) Lj (53) 

/=1 

where n = 1, 2 for T, L respectively and the integrals Jt,l, depending on /, are now calculated 
with 

i{m N rfAR A 



m 2 N V A QP[a(l - a) + m)/Q 2 fl 2 u [a(l - a) + m)/Q 2 f/ 2 
and the factor a(l — a) in (42) substituted as 

«(! " «) — M [a(1 r a) 2/n2 • ( 55 ) 
v ; a{l - a) +mj/Q 2 v ' 

The final EMC ratio is now calculated according to 

Evidently now it depends on two variables: zo defined by (54) and Q 2 . 

Calculating the proton structure function with (51) we find instead of (45) 

F 2p (x, Q 2 ) = {ll/f>\2)\N c i{Q/m N ) 2 -^ (57) 

where A is a factor resulting from the a and r integrations (and normalized to unity for 
(3 = 1 and mf = 0). To find it one has to calculate the integrals Jt,l with the function \{ z ) 
substituted by its lowest order term (l/3)z. If we denote the results J^\, then 



N 



f 



768 



ll7T 3 Z /=1 



£ + 2^ . (58) 
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It is evidently independent of z$ and can be calculated, say, with zq = 1. 
If we retain the form of which follows from (52) then (57) turns into 

q2-P(x) -A 

F 2p (x,Q") = c^ 7 =^ Y - , (59) 

with c a slowly varying function of lux, which we take a constant. Choosing at small but 
finite x the anomalous dimension in the form 

7 (x) = 2-/3(x) = l + a/lnx + b/ln 2 x, (60) 

we obtain a good agreement with all the experimental data on the proton structure function 
for x < 0.05 [14-17] with the parameters 

A = 0.237, c = 0.250, a = 3.87, 6 = 4.62 (61) 

(see Figs. 7,8). The formulas (59) and (60) then reasonably well extrapolate the asymp- 
totic pomeron to larger values of x, the power (3(x) taking into account the decrease of the 
anomalous dimension. 

Comparing (54) and (57) we observe that the nuclear structure function can again be 
written exclusively in terms of the proton one. Indeed the variable zq in (54) can be expressed 
similarly to (46) 

with an extra factor A in the denominator defined according to (58). 

Using (62) and (3 given by (60) and (61), we have calculated the EMC ratio for 
values of x < 0.05, Q 2 > 1 GeV 2 and A » 1, at which the experimental data are available 
[19-21]. Four quark flavours have been taken into account (the inclusion of the c quark results 
important). The value of the internucleon distance has been taken Rq = 1.35 fm. The best 
agreement with the data is achieved for values of the regulator mass m for light quarks of the 
order of 0.5 1.5 GeV. The results depend very little on m inside this interval: the overall 
change in 1 — does not exceed 10%. The results for m = 0.8 GeV together with the data 
are presented in the Table. As one observes, a very reasonable agreement is achieved, taking 
into account that the data, as a rule, are averaged over considerable intervals of x and Q 2 . 

At smaller x the behaviour of R^ with F2 P (x, Q 2 ) given by (59) naturally repeats the one 
obtained with (48) in the preceding section, except at low Q 2 where the effect of (effective) 
quark masses is felt. This is illustrated in Figs. 4-6, where the x-,Q- and ^4-dependence of 
R( A ) calculated with Fi v given by (59) is shown by dashed curves, to be compared with the 
one obtained with (48) (solid curves). 
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5 Conclusions 

As is well known there are serious problems of principle associated with the BFKL hard 
pomeron, related to the unitarity and the running coupling. In spite of all this, the hard 
pomeron approach remains the only one, based on first principles, which allows to study low 
x phenomena, although in a model theory. 

As we have found, the hard pomeron, at least, does not contradict the superlow x data for 
the proton structure function (x < 0.001). With some modification the pomeron description 
can be extended to the region of x < 0.05. This gives a motivation to apply it to the nuclear 
structure function at these x. 

With an additional assumption of small interaction between pomerons, which can be 
justified in the high N c limit, a closed formula is obtained for the nuclear structure function, 
expressing it directly in terms of the proton one. For a pure BFKL pomeron and a nucleus 
of constant density the EMC effect results a function of a single variable zq, which combines 
its x-, Q- and A-dependence. In principle zq can take on arbitrary values, so that can 
vary from unity to zero. However for realistic x, Q and A the variable zq is not so large, so 
that R( A ) remains not very different from unity, except for the lowest x and Q 2 and highest 
A. 

A novel feature of the hard pomeron approach is a changing anomalous dimension. As 
a result, the behaviour of the EMC effect on Q 2 depends on the value of x. For very small 
values of x the effect dies out as InQ/Q. However for larger x the Q-dependence becomes 
less pronounced until it practically disappears as the anomalous dimension becomes small. 

Comparison to the available experimental data for x < 0.05 shows a reasonable agree- 
ment. The agreement is not ideal, which may be explained partly by the imperfection of the 
existing data, averaged over considerable intervals of x and Q, partly by the presence of other 
subasymptotic terms different from the pomeron, which should be notable for not very small 
x. Evidently more data at lower x are needed to finally conclude about the relevance of the 
hard pomeron for the EMC effect. 
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7 Appendix: The asymptotics of the EMC effect at high Q 2 

With the quark mass mf retained and at high Q 2 the EMC ratio is a function of two 
small variables, = R( A \z ,l/Q 2 ). Since R^ A \d,l/Q 2 ) = 1 identically, the leading 

correction term is provided by R( A \zo, 0), that is, we can safely put rrif = to study it. 

The regions of small a and I — a are responsible for the leading behaviour at zq — > 0. 
Therefore we can limit ourselves with the transversal cross section and consequently with the 
integral Jt, Eq. (36). Transforming to the integration variable w by 

a(l-a) = l/{w 2 + 4), (63) 

the dominant contribution will come from large w. Choosing wq » 2, we find that at small 

roo roo 

J T ~ 2 / rdrK?(r) / — X (*), (64) 

where 

z = z (wrf, 1<0<2. (65) 
Passing to the integration over z, we present the right-hand side as 

, roo roo A 7 

J T ~(2/(3)z 2 /(S j Q r*drK 2 (r)j^ X (z), (66) 

with z\ = zoiwor) 13 . 

The internal integral over z can be transformed into 



roo 



dz 



-2/(3 



where 



/oo 
dxx~ 2//3 exp(-zx) (68) 

can be expressed via the incomplete gamma function. Evidently the first term in (67) gives a 
contribution to Jt regular st zq = and so a correction to R^ of the order zq ~ 1/Q. The 
leading correction comes from the term with F{z\). Straightforward estimates give that at 
small z 

F(z) ~ a + bz 2 ^' 1 , l</3<2; 

F(z) ~ a + bzlnz, 0=1. (69) 
Putting this into (67) and then into (66) we obtain our final result. At high Q 2 

R^ ~ 1 - c Q 1_2//3 , 1< < 2; 

#(^) ~ 1-chxQ/Q, 0=1. (70) 

In particular, with a dimension close to its canonical value = 2 the EMC effect results 
independent of Q. 
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Table 

EMC ratios at experimentally known points. 



A 


X 


Q 2 GeV 2 


R 


Rexp 


Ref 


131 


0.0065 


1.34 


0.780 


0.840± 0.107 


[19] 


131 


0.0095 


1.91 


0.825 


0.800± 0.108 


[19] 


131 


0.0210 


3.58 


0.889 


0.900± 0.110 


[19] 


64 


0.0310 


4.50 


0.925 


0.930± 0.033 


[20] 


64 


0.0500 


8.50 


0.955 


0.955± 0.028 


[20] 


64 


0.0150 


4.50 


0.923 


0.857± 0.032 


[20] 


64 


0.0310 


3.30 


0.904 


0.963± 0.022 


[20] 


64 


0.0500 


6.40 


0.944 


1.005± 0.018 


[20] 


12 


0.0055 


1.10 


0.861 


0.904± 0.013 


[21] 


12 


0.0085 


1.60 


0.894 


0.939± 0.012 


[21] 


12 


0.0125 


2.20 


0.917 


0.939± 0.009 


[21] 


12 


0.0175 


2.90 


0.934 


0.957± 0.010 


[21] 


12 


0.0250 


3.60 


0.945 


0.963± 0.009 


[21] 


12 


0.0350 


4.50 


0.954 


0.990± 0.009 


[21] 


12 


0.0450 


5.50 


0.962 


0.983± 0.010 


[21] 


40 


0.0085 


1.40 


0.840 


0.846± 0.014 


[21] 


40 


0.0125 


1.90 


0.871 


0.870± 0.011 


[21] 


40 


0.0175 


2.50 


0.895 


0.908± 0.011 


[21] 


40 


0.0250 


3.40 


0.917 


0.946± 0.009 


[21] 


40 


0.0350 


4.70 


0.937 


0.956± 0.009 


[21] 


40 


0.0450 


5.70 


0.946 


0.986± 0.011 


[21] 
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Figure captions 

Fig. 1. Pomeron diagrams for the EMC effect in the leading approximation in l/N c . With 
a small coupling constant the diagram (a) dominates at large energies. 
Fig. 2. EMC effect as a function of the parameter zq (Eq. (44)). 

Fig. 3. The pomeron fit (48), (50) to the proton structure function for x < 0.001 compared 
to the experimental data [14-17] at Q 2 = 3, 12 and 35 GeV 2 . 

Fig. 4. The x-dependence of the EMC effect calculated with the proton structure function 
described by the pure pomeron (Eq. (48), solid curves) and the subasymptotic pomeron (Eq. 
(59), dashed curves) for x < 0.001, A = 64 and Q 2 = 3 and 50 GeV 2 . 

Fig. 5. The Q 2 -dependence of the EMC effect calculated with the proton structure function 
described by the pure pomeron (Eq. (48), solid curves) and the subasymptotic pomeron (Eq. 
(59), dashed curves) for A = 64 and x = 10 -3 and 10~ 6 . 

Fig. 6. The yl-dependence of the EMC effect calculated with the proton structure function 
described by the pure pomeron (Eq. (48), solid curves) and the subasymptotic pomeron (Eq. 
(59), dashed curves) for x = 10" 3 and 10~ 6 and Q 2 = 10 GeV 2 . 

Fig. 7. The subasymptotic pomeron fit (59)-(61) to the proton structure function for 

x < 0.05 compared to the experimental data [14-17] at Q 2 = 3, 12 and 35 GeV 2 . 

Fig. 8. The effective anomalous dimension 7 = 2 — (3 as a function of x (Eqs. (60), (61)). 
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